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Two particle correlations are studied in the reaction plane ol peripheral relativistic heavy ion 
reactions where the initial state has substantial angular momentum. The earlier predicted rotation 
effect and Kelvin Helmholtz Instability, leads to space-time momentum correlations among the 
emitted particles. A specific combination of two particle correlation measurements is proposed, 
which can sensitively detect the rotation of the emitting system. Here the method is presented 
in simple few source models where the symmetries and the possibilities of the detection can be 
demonstrated in a transparent way. 
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I. INTRODUCTION 

Collective flow is one of the most dominant observable 
features in heavy ion reactions up to the highest avail- 
able energies, and its global symmetries as well as its 
fluctuations are extensively studied. Especially at the 
highest energies for peripheral reaction the angular mo- 
mentum of the initial state is substantial, which leads 
to observable rotation according to fluid dynamical es- 
timates [1 . Furthermore the low viscosity quark-gluon 
fluid may lead to to initial turbulent instabilities, like 
the Kelvin Helmholtz Instability (KHI), according to nu- 
merical fluid dynamical estimates [2J, which is also con- 
firmed in a simplified analytic model [3] . These turbulent 
phenomena further increase the rotation of the system, 
which also leads to a large vorticity and circulation of the 
participant zone. [3]. In ref. [5] it is estimated that the 
increased rotation can be observable via the increased v±- 
flow, but the i>i signal at high energies is weak, so other 
observables of the rotation are also needed. 

The two particle correlation method is used to deter- 
mine the space-time size of the system emitting the ob- 
served particles, thus providing valuable information on 
the exploding and expanding system at the freeze out 
stage of a heavy ion collision. This method is based on 
the Hanbury Brown and Twiss (HBT) method, originally 
used for the determination of the size of distant stars [5] . 
In heavy ion collisions the HBT method was used first for 
the same purpose, the determination of the system size 
[S] , but later also the ellipsoidal shape of the system and 
its tilt [5] ■ It was also observed relatively early that 
the expansion of the system modifies the size estimates 
due to the collective radial flow velocity of the emitting 
system [S], while the effect of flow on two particle corre- 
lations was also analysed in great detail [ID]- Transport 
model studies have indicated that the HBT radius shows 
a minimum at the phase transition threshold [llj . 

A detailed two particle correlation study of flow ro- 
tation was not performed up to now [12] . This became 
actual as at higher beam energies, where the initial an- 
gular momentum of the participant system is increasing 
in peripheral reactions, the system may rotate, causing 
a significant and detectable effect. Here we do not want 



to use the most advanced state of the art developments 
of the HBT method, rather we want to show in the most 
simple way, how rotation can be detected by the method. 

Different theoretical approaches were worked out up to 
now to evaluate the two particle correlation functions in 
different reaction models. We here picked one method, 
which is not used very frequently up to now, but it can 
be generalized and used to fluid dynamical models with 
arbitrary flow patterns well. With this method hereby we 
study simplified, idealized fluid dynamical systems with 
different symmetry structures. These studies show how 
we can detect rotation via two particle correlation func- 
tions and what effects may cause difficulties in identifying 
rotation. 

Based on these studies we also present a Differential 
Hanbury Brown and Twiss (DHBT) method, which can 
sensitively determine the strength (vorticity or circula- 
tion) of the rotating flow and the direction of this rota- 
tion. 



A. The Emission Function 

Following [T3] the definition of a particle 4-current is: 

N" = J P^f(x, P ) , (1) 

where f(x,p) is the invariant scalar phase space density 
distribution of the emitted particles. The total flow of N 
nucleons across a space-time (ST) hypersurface, with the 
surface element, da^ is: 

N = J d ^Jda,p»f(x,p). (2) 

If we do not perform the integration over the momentum, 
p, then we get the Cooper-Frye formula |14j for invariant 
momentum distribution: 

E~ = J da^f(x,p) , (3) 
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where E — p°. This assumes that there is a 3-dimensional 
Freeze Out (FO) hypersurface, which can also be gen- 
eralized to become a layer, as we will show. Using 
d 3 p = dp\\dpj_ = p°dydp^) we can write this in another 
form, and we can extend it to a 4-volume integral of a 
4-dimensional Source Function, S, as [T5j [16] : 



dN 



dyd 2 p A 



da^(x)p l± f(x,p) = J d 4 xS{x,p) 



(4) 



where we assumed that the emission appears on a 3- 
dimensional hypersurface with the outward pointing nor- 
mal, <i<7 M . The emission function gives the distribution of 
the ST positions of momenta of emitted particles. The 
emission function gives the number of the particles, AN, 
emitted in the phase-space element A 3 xA 3 p per unit 
time, At. We get a Lorentz invariant scalar if we also 
multiply it by the energy, p°, of the emitted particles: 



S(x,p)=p° 



AN 
AtA 3 xA 3 p ' 



(5) 



Along the lines of this introduction we can describe the 
emission over a hypersurface also as a 4-volume integral. 
We still assume a ST hypersurface with an outward point- 
ing normal vector, dcr M , but then the ST integral is in- 
terpreted as a 4-volume integral with a delta function for 
the given surface, 5{x' — x) as: 

J S(x,p)d 4 x = J d 4 xp^a t ,(x')5^(x' -x)f(x,p) , 

(6) 

where the emission is constrained to a 3-dimcnsional hy- 
persurface in the ST and it is directed in a given direction 
characterized by the unit vector in the source described 
by &ij,(x') S^'(x'). Even if we assume that the source is 
in a ST layer (which is not too thick, e.g. 2-3 fm or 
2-3 fm/c), in this layer we can have a maximum of the 
emission. This FO layer in case of pions is narrow if wc 
consider the rapid and simultaneous hadronization and 
FO from the plasma. This could even be idealized as a 3- 
dim hypersurface in the 4-dim ST if the thickness of the 
layer is neglected. Thus, <7 M , is the unit normal vector of 
this surface or layer can be both timelike or spacelike: 

(T^fTp = +1 for timclikc hypersurface or layer , 
cr^cTp = —1 for spacelike hypersurface or layer . (7) 

The idealization of FO in a 3D hypersurface is not 
necessary, however it makes the presentation more trans- 
parent. The more realistic emission distribution must 
happen in a 4D ST surface layer, which can still can 
have an effective space-like or time-like normal vector. 
See refs. [HI ] IT7HI9] . Although, one might naively believe 
that in case of a normal vector, = (1, 0, 0, 0) the emis- 
sion is uniform in all spatial directions, this is not true, as 
the local flow velocity also influences the emission prob- 
ability [101 lIHj- When the flow velocity points in the 
direction of the detector (or the FO normal) the prob- 
ability for the emission into the a^-direction is bigger, 



so that the emission probability should be proportional 
to (w M • cTfj). This will be important later on when the 
observability of rotation is discussed. 

At ultrarelativistic energies most of the FO happens in 
time like directions, frequently idealized as t ^constant 
or r =constant hypersurfaces or layers. In addition even 
in case of time-like FO the deeper (or earlier) points of 
the FO layer have a smaller emission probability because 
of the opacity of QGP, indicated also by the strong jet 
quenching. This effect causes additional asymmetries in 
the emission. 



B. Hydrodynamical Parameterization 

Let us assume that at the points of the source the 
matter is still in local equilibrium. Then we can describe 
the phase-space distribution of the particles, f(x,p) by a 
Jiittner distributions [13) or a relativistic Bose-Einstein 
or Fermi-Dirac distribution. Furthermore, instead of the 
delta function we may assume the emission distributed 
in a ST layer, which still has a preferred direction of 
emission . 

For our purposes the most suitable parametrization of 
the emission function is introduced in the special " Buda- 
Lund" model, see section 8 of ref. [20J . Here the emission 
function is parametrized as 



S(x,p)d 4 x=p»d 4 Z f ,(x)f(x,p), 



(8) 



where the emission probability and its dependence of the 
FO direction is already included in the p^ d 4 Yj^{x) term. 
The 4-volume integral is directed and yields a maximum 
for which is closest to d 4 Y,^(x). In the Buda-Lund 
model the FO direction points into the flow 4-velocity, 
d 4 T, fl (x) cx vf 1 . This is also a frequent approximation in 
other fluid dynamical models, although, in the general 
case it is not a valid approximation. Thus, this approach 
is identical to the dynamical volume FO in a layer |17l - 
IT9] . discussed above. We also assume that the FO layer 
depends parametrically on the proper time, r, (or dis- 
tance s) in the FO factor d 4 E^(x) oc a fi d 4 x. So, that the 
emission probability is proportional to G(x) H(t): 

S(x,p)d 4 x=p fl <5>(x) G{x)H(t) drd 3 x f(x,p) , 



where 



H(r) 



(2 7 r9 2 ) 1 /2 



exp 



(r-f) 
29 2 



2 1 



(9) 



and that the widths of the emitting sources do not change 
significantly during the course of emission from a given 
source. (The emission probability will be discussed in 
more detail when the source asymmetry is analysed in 
section [TTT] ) The emission function characterized with a 
locally thermalized volume-emitting source is then: 



S(x,p) d 4 x 



g p> 1 a fl (x)G(x)H(T)dTd 3 x 



exp 



V ^(x) 



T(x) 



(10) 
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where in place of f(x,p) we inserted the relativistic Bose- 
Einstein distribution. The factor g stands for the de- 
generacy, u^(x) is the 4- velocity field, T(x) is the tem- 
perature field, n(x) is the chemical potential and d^x = 
dr dx dy dz. G(x) is the ST emission density across the 
layer of the particles (e.g. pions). This can be approxi- 
mated with a ST hypersurface and then with — > we 
obtain the Cooper-Frye FO description. [T4"] 

For the phase space distribution we frequently use the 
Jiittner (relativistic Boltzmann) distribution: 

which is normalized to the invariant scalar density of par- 
ticles 

n(x) = N»u» =u,j f(x,p) = ^C n , (12) 

where C n — Anra 2 TK2{rn/T) and we use the c = k = 1 
convention. Thus in terms of the local invariant scalar 
particle density the Jiittner distribution is |13j 



j n(x) ( p^u^x) 



T(x) 



(13) 



We can also define do* 1 = dx dy dz a* 1 (for a timelike 
surface or layer) where the norm of da^ is the 3-volume 
of the source element (like a fluid cell), similarly to Refs. 

TOM 



Our source function in this case, with Eq. (13), sim- 
ilarly to Ref. [ID] in the frame where <r M = (TO, 0, 0) 
is: 



G(x) = N FD (x)/n(x), 



(14) 



where Nfd (x) is the density of particles arising from 
fluid-dynamical (FD) evolution (or from other transport 
models). Then taking inti account the flow velocity in 
this frame: 



Nfd(x) = 771(35) 



(15) 



The H(t) (or H(s)) freeze-out probability along the r 
(or s) parameter, across the layer can be integrated sep- 
arately from the remaining three orthogonal coordinates. 

Here we assume that the primary direction of the emis- 
sion is <T M . Thus, in case of an explosively expanding 
system it points towards the detector, so k^. The 

emission happens from the 4-volume of ST surface layer 
with an effective normal direction er^ 1 , and not from a ST 
hypersurface. 

In ref. |10j the correlation function was analyzed in 
detail in dependence of the direction of the primary emis- 
sion direction, fc M . A possibility of longitudinal momen- 
tum difference was considered in terms of the rapidity of 
the two emitted particles, where the difference of the lon- 
gitudinal momenta, the width parameter, was studied in 
detail. While this analysis provides a deeper insight into 



the features of the correlation function, in our studies 
we concentrate to the rotation of the system and restrict 
oursclf to a simplest presentation of the correlation func- 
tion, which can be realized experimentally without much 
additional effort. 

We can assume FO from a narrow layer at a proper 
time hyperbola Tpo — cons, like in the Buda-Lund 
model, see ref. [20J. This can be practical if the CFD 
model uses proper time and rapidity coordinates. Re- 
cent studies indicate that irrespective of the coordinate 
system choice, in the major part of FO in high energy 
collisions the FO happens near to a constant proper time 
hyperbola, although the origin of this FO-hyperbola is 
at an earlier point of time than the intersection of the 
centers of the projectile and target trajectories [2Tj . 

For the first test purpose we take an oversimplified 
model of 4 fluid elements, which may or may not expand 
or rotate. We assume that these are in the reaction plane, 
[x — z]-plane, and will characterise parameters of a heavy 
ion reaction based on CFD results. We assume that the 
system is stationary so the time emission probability is a 
Gaussian (like) distribution in time. Later on we intend 
to study and see that these methods can be applied for 
realistic full scale FD calculations also. 



C. Pion Correlation Functions 

The pion correlation function is defined as the inclusive 
two-particle distribution divided by the product of the 
inclusive one-particle distributions, such that |15j : 



C(pi,P2) 



Pl( P l)Pl(p2Y 



(16) 



where p\ and P2 are the 4-momenta of the pions. 

We will assume pions are created at two points, X\ and 
X2, which are distributed in space. The particle distribu- 
tion is given by the reduced phase space source distribu- 
tion: 



S(x,p) = f(x,p)pf l a ft G(x)H(T) 



(17) 



For two identical pions with momenta p\ and P2 the two- 
particle distribution is: 

P2(Pi,p2) = / d 4 x 1 d 4 x 2 S(xi,pi)S(x2,P2)\i>i2\ 2 , (18) 



(19) 



where the wave equation i/>i2 is given by: 

' ' ip 1 -x 1 +ip 2 -x 2 _|_ e ipi-x 2 +ip2-X!^ 



We now introduce the center-of-mass momentum 



k = +P2), 



(20) 



1 The vector k is the wavenumber vector, k = p/K so for numerical 
calculations we have to use that he = 197.327 MeV fm. 
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and the relative momentum 

Q=Pi ~P2, 



(21) 



where assuming the mass-shell constraint for the two par- 
ticles and so we have q-k — (pi — p 2 ) ■ (pi +P2)/2 = (pi — 
pl)/2 = (ml - ml)/2 = 0, which leads to q° = q ■ k/k° 
With the relative and center-of-mass momentum we 
can write the wave equation as: 



ik-(x 1 +x 2 ) 



12 



V2 



e iq-{x 1 -x 2 )/2 _|_ e -iq-(xi-x 2 )/2 



and then 

2 



_|_ 1 ^gig-(xi-x 2 ) _|_ e -iq-(xi-x 2 ) S j 



(22) 



(23) 



We can then insert Eq. (23)intoEq. ( 18 ), and we obtain 



ft(pi,Pa)= / d i x 1 d i x 2 S(x 1 ,k + q/2)S(x 2 ,k-q/2) 



> (24) 



where the last term in the brackets is cos[g(a;i — x 2 ) 
Similarly for the one particle distribution we get: 



P x (p) = / d A x S(x,k) . 



(25) 



We now use a method for moving sources presented in 
ref. [32]. Using Eqs. ( 24|25 ), together with the definition 
of the correlation function we have: 



C(k,q) 



R(k,q) 



|/ d 4 xS(x,k)\ 



(26) 



where 



R(k,q)— J d 4 xi d 4 x 2 coa[q(xi — x 2 )]'x 
S(x 1 ,k + q/2)S(x 2 ,k-q/2) . 
Here R(k,q) can be calculated [52] via the function 

J(k,q) — J d^x S(x,k + q/2) exp(iqx) = 
'd^S( X ,k + q/2)[ C os(qx) + iM^], 



(27) 



(28) 



and we obtain the R(k,q) function as 

R(k,q) =Re [J(k, q) J(k, -q)) 



(29) 



This can be verified, by using Eq. ( 28 ), forming a dou- 



ble integral over d A x\ d 4 x 2 from J(k, qjj (k, —q), yielding 
to a term exp[— iq(x\ — x 2 )]. Then taking the real part 
of the double integral leads to a term cos[g(xi — x 2 )\ and 
this recovers Eq. (27). 



Source with local Jiitner distribution: Let us take 



the S(xi,pi)S(x 2 ,p 2 ) term in Eq. (24), and assume that 
the single particle distributions, f(x,p), in the source 
functions are Jiittner distributions, which depend on the 
local velocity, u^(x) , via the term 



exp 



-p^u^x) 
T(x) 



(30) 



as shown in Eq. ( 13 ). Here the local flow velocity may be 



different in different locations, x\ and x 2 , and this influ- 
ences the correlations of the observed momenta. Thus, 
the scalar products in terms of k and q become: 

exp(-piiti) exp(-p 2 "2) = 

exp(-(k + q/2) ui) exp(-(fc - q/2) u 2 ) = (31) 
exp(— ku\) exp(— ku2) exp(— q(u\ — U2)/2) 

where we used the notation u\ — u(x x ) = ^(xi). We 
assume that for a given detector position the normal di- 
rection of the emission is approximately the same, so for 
the two sources the term p^a^ (x) is the same and it can- 
cels in the nominator and denominator. 

Thus, the expression of the correlation function, Eq. 



(271 will be modified to 



R(k,q) — d 4 xid 4: x 2 S(xi, k)S(x 2 , k) cos[q(xi—x 2 )}x 



exp 



u(xij u(x 2 ) 
T(x t ) ~ T(x 2 ) 



(32) 



and the corresponding J(k, q) function will become 

q ■ u(x) 



J(k,q) = d x S(x,k) exp 



2T(x) 



cxp(iqx) 



In Eq. ( 32 ) the term including the momentum compo- 



nent q and the flow velocity u becomes unity if the source 
has a uniform distribution of u(x)/T(x), and in this case 
we may be able to use the so called smoothness approx- 
imation: S(x,k+q/2)S(y,k—q/2) S(x,k)S(y,k), and 
the correlation function, expression ( 26 ) , takes the form 



C(k,q) = l + 



|J d A x e iqx S(x,k)\ 
\J d 4 xS(x,k)\ 2 



(34) 



D. One Fluid Cell as Source 

We now assume a source function, which is reduced to 
one Freeze Out (FO) time moment. Thus the integra- 
tion over the 4-volume of an emission layer is reduced 
to the 3-volume of a FO hypersurface. For simplicity, 
we assume FO along the timelike cordinate, t, where we 
assume a local Jiittner distribution. Thus, we have the 
source function as 



,V (l ./,) = (V(,!//(/!, M )(-^^ 



fc M <J M , (35) 
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where k^a^ is an invariant scalar, and for a single cell we 
use a simple quadratic parametrization for n(x) as: 



G(x) = jn(x) = 7n s exp — 



y 2 + z 2 



2R 2 



(36) 



Here n s is the average density of the Gaussian source (or 
fluid cell) of mean radius R. 

Single source at rest: The invariant scalar k^u^ can 
be calculated in the frame where the cell is at rest. We 
have then 

^ = (1,0,0,0)^-^ = -^ = -^. (37) 

In this simplest case we also assume that the FO direction 
is = (1, 0, 0, 0), so the r-coordinate coincides with the 
t-coordinate, and it is orthogonal to the x,y,x— coordi- 
nates. Then we can make use of the following integral: 



'd A x = 



(38) 



We can perform the integral along the t direction of H(t), 
which gives unity and then the single particle distribution 
is 



d 4 x S{x,k) = ^ (fc M (7 M ) exp \ ~^r ' 



+oo r+oa ^ 2 r+oo ^ 2 

H{t)dt I e~?R I dx I e~2R^dy 



e 

-oo 



z 2 

2n7 dz = 



n s (k^a^) exp 



2irR 2 



(39) 



where T s is the temperature of the source, and Ef. = fc° 
in the rest frame of the fluid cell. Due to the normaliza- 
tion of H(t) the integral over the time t is unity. The 
contribution to the nominator from Eq. (|33|) is 



J(k,q)= I d*xe i **e-'flV r ''>S{x,k) = n ' ( f^- x 



exp 



E k +q°/2 



+ 00 2 

— y 

e 2H 2 e 

— oo 



n s {k^a^) 



iq y y dy 

>2N3/2 



+oo ^ 2 

e^fii 1 e~ lq:cX dx x 

— oo 

+ 00 ^ 2 

e^fa 1 e~ lqzZ dz = 



H(t)e lq0t dt 



exp 













exp 


2T S _ 


X 





R 2 n 




exp 


2 q 


exp 



e 2 

2 



where we used /_ exp(— p 2 x 2 ±qx)d. 



(40) 



exp(q 2 /(4p 2 )) t 2 3 3.323/2. In the time integral the 
present choice of would give (q ) 2 , but we wanted 



to indicate that other choices are also possible and they 
would yield (a^q^) 2 . In the J(k,q)J(k,—q) product the 
terms exp[±<7°/(2T s )] cancel each other. Inserting these 
equations into ( 26 ) we get 



C{k, q) = l + exp (-e 2 (<7%) 2 ~ R 2 1 2 



(41) 



If we have a source at a point in the FO layer, which 
is at a longer distance from the external side of the FO 
layer than O, then the contribution of the time integral 
from this point is reduced. In a few source model it is 
more transparent to describe this reduction by assigning 
a smaller weight factor to the contribution of the deeper 
lying source. See in section [Tn] 

If we tend to an infinitely narrow FO layer, O — > 0, i.e. 
to a FO hypersurface, then 



C(k, q) = l + exp (-R 2 q 2 ^ 



(42) 



The k dependence drops out from the correlation func- 
tion, C(fc, q) as the k dependent parts are separable. See 
Fig. [l] The size of the fluid cells in a high resolution 
3+1D fluid dynamical calculation is (0.3fm) 3 . With this 
resolution the numerical viscosity of the fluid dynamical 
calculation [24 is the same as the estimated minimal vis- 
cosity of the QGP [5S] which occurs at the critical point 
of the phase transition [26]. As Fig. [I] shows the correla- 
tion for such a cell size yields to an extended distribution 
in the relative momentum q. 




1 2 

q (1/fm) 

FIG. 1. (color online) The correlation function, C(k,q), for 
a single, static, spherically symmetric, Gaussian source with 
different radii, R = 4, 1 and 0.25 fm, (blue dotted, red dashed, 
and full black lines respectively), as described by Eq. (42 l. 



For the study of the rotation of the system the thick- 
ness of the FO layer is of secondary importance, espe- 
cially if we discuss only a few fluid sources. In this 
case the role of the depth of a source point within the 
layer is given by its reduced contribution to the particle 
emission. This can be represented much simpler with as- 
signing emission weights to the small number of sources. 
Thus, in the following discussion, we do not go into the 
details of the time structure of the emission. 
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Single moving source: Let us take a single source 
which moves in the x-direction with a velocity v x . Then 



we have, 



7s (1 7 v x , 0, 0), and the scalar product 



k ■ u s /T s = k^u^/Ts provides an additional contribu- 
tion to the correlation function. However, in the case 
of a single fluid cell or a single source the velocity and 
the temperature do not change within the cell, so the 
modifying term in eq. ( 32 ) becomes unity. We use 
j(Ek — k x v x ), and the source function becomes 



k^ 



S(x,k) = ^ exp 

<->rt 



k ■ U s 



(43) 



where [k ■ u s /T s ] = ["/ s (E k - k x v x )/T s }. 

Within the source (or fluid element) the velocity u s 
and temperature T s are assumed to be the same. The 
source or fluid clement may have a density profile, but 
this profile should be the same for all cells (although the 
average density, n s is not the same for all cells. The 
spatial integrals can be performed in the rest frame of 
the cell, giving the same integral result as above (39 1, 



because the moving cell-size shrinks, but the apparent 
density increases, so that the total number of particles in 
a cell remains the same as it is an invariant scalar. 
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FIG. 2. (color online) Two steady sources in the reaction 
([a; — z]) plane with a distance between them of 2d in the 
a;— direction. 

Then the integral of the single particle contribution is 

d 3 xS{x,k) = 



n s {k^a^) exp 



k ■ it. 



2nR 2 



(44) 



Then the two particle distribution: 

J(k, q) = J d 3 x e lq - x S{x 7 k) exp 
k ■ u. 



n^fc^oyjexp 



(2irR 



T s 

2\ 3 / 2 



exp 



exp 



q ■ u s 
q ■ u s 

2T S 
R 2 9 



(45) 



When calculating R(k, q), in the J(k, q)J(k, —q) product 
the terms exp[±q ■ u s /(2T s )] cancel each other. In the 



formulae the h = 1 convention is used and k and q are 
considered as the wavenumber vectors. 



We then insert these equations into equation ( 26 ) and 
we get for one moving Gaussian source 



C(k,q) = 1 + exp (-R 2 q 2 ^ 



(46) 



Again, this result does not depend on k, just as the pre- 
vious single source at rest, Eq. (42). 



II. SYMMETRIC FEW SOURCE MODELS 

A. Two Steady Fluid Cells 

For emission from two steady sources, two particle cor- 
relations were studied in ref. |20j . Here we use the 
present method. We assume that the two source system 
is symmetric both their positions are placed symmetri- 
cally and also their FO normal vectors, <r M , are the same. 
If the normal were <r M = (1,0,0,0), then the invariant 
scalar fc M <7 M would be k° — Ej~, although we do not need 
this additional requirement to illustrate the correlation 
function, which would arise from an idealized symmetric 
system. 

We also assume that the time distributions, H{t) for 
the two sources are identical, so these can be integrated 
simultaneously and yield unity. If we have two sources 
then the source function is 



S(x,k) = }^S s (x, k) 

S 



k • u s 



(47) 



while the J function in the Jiittner approximation is 
J(k,q) = 



H eX P 



q ■ u s 



2Ts 



exjp(iqx s ) I d x S s (x,k) exp(iqx) , 
's 

(48) 



where x s is the position of the center of the source, and 
the spatial integrals run separately for each of the iden- 
tical sources, i.e. we assume fluid cells with identical 
density profiles, but with different densities, n s and tem- 
peratures, T s . 

In case of steady sources = (1, 0, 0, 0), and the spa- 
tial integral for one source is the same as for a single 
source. Thus, 



J d 3 x S(x, k) = ^2j d 3 x S s (x, k) 



2nR 2 f 2 (fc^)^-^exp(- 



Ek 



(49) 
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and 
J(k,q) = 
^exp 



2T, 



exp(iqx s ) / d x S s (x, k) exp(iqx) 
's 



exp(iq°x®) exp(—iqx s ) 



exp 



exp 



g 

2T< 



(50) 



In the J(k,q)J(k,—q) product the terms exp[±g°/(2T s )] 
cancel each other. Both J(k, q) and J(k, —q) include a 
sum [exp(iqx s ) +exp(—iqx s )], and their product leads to 
a factor 2[l+cos(2gcc s )]. Here we assumed that the time- 
like extent of the emission layer is negligible compared to 
the space-like size. 

Consequently, if the two sources have the same param- 
eters, just different locations, x\ = — X2 (see Fig. [2]) then 



C(k,q) = 1+ -cxp(-i?V)[l + cos(2ga; s )] (51) 




1 2 

q (1/fm) 



FIG. 3. (color online) The correlation function, C(k,q), for 
two spherical, Gaussian sources at rest, shown in the direc- 
tion of the displacement (here q x ) with different distances, 
d — 2.5, 1.25 and fm, (dotted blue, dashed red, and full-black 
lines respectively) and in the orthogonal direction, q y which 
is identical with the d = contribution (full black line). For 
finite d the cosine term results in a modification in the direc- 
tion of the line joining the sources, which becomes apparent 
if the distance between the sources becomes small compared 
to the size of the spherical sources. 



only. In the other directions, q y and q z , the single source 
result (42) is returned. 



If the distance of the two sources, 2d, is comparable 
or smaller than the radius of a single source, R, then 
the two source configuration leads to visible zero points, 
C(k, q) — 0, on the q x -&xis at 2q x d = ±(1 + 2n)ir, where 
n = 0,1,2,3,... . In Fig. [3] for the d = 2.5fm case we 
see these zero points at q x — n/(2d), 2>ir/(2d) , ... , while 
at the points q x = 2ir / (2d) , Air / (2d) , ... the distribution 
function, C(k,q x ) touches (becomes tangent to) the dis- 
tribution function for d = or the distribution function 
C(k,q y ). 

The appearance of the zero points is to a large extent 
an artifact of the used very simplistic two source model. 
In case of other additional sources these zero points would 
disappear. Nevertheless, this feature illustrates that the 
correlation function can be more complex than a set of 
Gaussians of the momentum difference q in different di- 
rections or at different rapidities. 



t- 



FIG. 4. (color online) Two moving sources in the reaction 
([x — z\) plane with a distance between them of 2d in the 
a;— direction. The sources are moving in the directions indi- 
cated by the (red) arrows. 



B. Two moving sources 

We study the system the same way as before, but now 
the two sources are moving in opposite directions, so 



that u s = tii or t/2 where 
(isHsi-V!)), and u s = 
Fig. HI Similarly, x s = 
x% = (t s , —x s ), and x\ = 
that FO happens at a t 



= (ls,7sVl), u 2 = < = 

7 S v si so that Ui = —U2, see 
X\ or xi where x^ = (t s ,x s ), 
—x 2- For now we also assume 
=const. FO hypersurface, so 



da^ = (1, 0, 0, 0) and so h=t 2 . 

If we have several sources then the source function in 
Jiittner approximation is 



This result agrees with ref. [2D]j section 9.1 (p. 41), 
and in the limit of x s = it returns the single source 
result, Eq. (42 1. See Fig. [3] If the distance of the two 
sources is 2d, i.e. xi = d and x 2 — —d, then 2qx s = 
2q x d, thus the modification appears in the g^-direction 



S(x, k) = ^2 S s (x, k) = 

s 

,, ,, . . Tl 3 (x) 

(^m) Ett 1 exp 



k- 



(52) 



while the J function is 



J(k,q) =J^exp 



q ■ u s 

2T S 



(53) 



exp(iqx s ) / d x S s (x,k) exp(iqx) , 
's 



where x s is the 4-position of the center of source s, 



and the spatial integrals run separately for each of the 
identical sources, i.e. we assume fluid cells with identical 
density profiles, but with different densities, n s , veloci- 
ties, u s and temperatures, T s . 



The spatial integral for one source is the same as for a 
single source. Thus, 



J d 3 x S(x, k)=J2j d * x S.(x, k) = (Jfc" (2ttR 2 



3/2 n s 

^ n 



exp 



exp 



ku s 



exp 



ku 6 



(54) 



This returns Eq. (49) if = (1,0,0,0). The function J(k,q) becomes 

q ■ u 



271 



J( k ,q) = E exp 

S 

(fc M <r M ) (27ri? 2 ) 3/2 exp 



exp(iqx s ) j d x S s (x,k) exp(iqx) 
s 



R 2 q 2 \ 
2 J 


£- 

s 


R 2 2 \ 


n s 







n s 


k ■ u s 




q-u s ~ 


7^ ex P 


Ts 


exp 


2T S 





ku s 




qu s ' 


exp 


[ T s _ 


exp 


2T S 



exp(— iqx s ) + exp 







r q° i 
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T s 


exp 


2 T s _ 


ku s 




qu s ~ 






exp 


2T S _ 


e 



exp(iqx s ) 

cxp(iq°x° s )x 
cxp(iqx s ) , 



(55) 



where the factor exp(iq°x®) ca n be dropped if the FO time distribution is simultaneous for the two sources, because 
then x° s = 0. This returns Eq. d50) if < = (1, 0, 0, 0). 



Now we can divide the two particle correlation with the square of the single particle distribution 

2ku, iku. 



Re [J(k,q)J(k,-q)] 
\f d A xS{ Xl k)\ 2 



Re 



exp(-i?V 



qu s qu s 

e ^r^ e i2qx s +e -*Tr- e -i2qx, 



ku, ku, 

e T » + e T = 



2„2\ 



exp(— R q 



Re 


2cosh( 2 ^ M ^ 


qu s qu s 

+ e T ° (cos(2qx s ) + isin(2qx s )) + e T » (cos(— 2qx s ) + i sin(— 2qx s )) 


2 




+ 1 





exp(-R 2 q 2 



cosh + cosh (%gA cos(2qx s ) 



cosh 



/ 2ku 3 \ 
V T * ) 



(56) 



Consequently, if the two sources have the same param- 
eters, just opposite locations with respect to the center, 
and opposite velocities, then the correlation function is 



if q = 0. 



C{k,q) = 1 + exp(-R 2 q 2 ) 



cosh (^fA + cosh (^A cos(2qx s ) ( 57 ) 
cosh (^) 

This returns Eq. ^ if < = (1, 0, 0, 0), and C(k, q) = 2 



If we have two sources placed at x — ±d x , and with 
the velocity in the indirection, ±v z , then the correlation 
function is for the different directions becomes: 
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C(k x , q x ) =1+2 exp(-i? 2 ^) [1 + cos {2q x d x )\ 
C(k x ,q y ) = 1 + exp(-I? 2 g 2 ) 



C(k x ,q z ) = 1 + 



exp(— R 2 q 



1 + cosh 



(58) 



C(fc y , & ) = 1 + - exp(-i? 2 gr 2 ) [1 + cos (2^4,)] 
C(k y ,q y ) = 1 + exp(-i? 2 g 2 ) , 

C(/c< / ,^) = l + ^ — 1 + coshl^- 



(59) 




1 2 3 

q (1/fm) 

FIG. 5. (color online) The correlation functions, C(k,q), for 
two moving sources where the displacement of the sources is 
in the ^-direction, and the center-of-mass momentum, k, of 
emitted particles is in the x and y direction. The solid black 
line is for the momentum difference, q z , the dashed red line is 
for q y and dotted blue line is for q x . The radius of the sources 
is R — 1 fm, (same as in Fig. j3j), the displacement is d = lfm, 
and the source velocity is, r yv z /T a = 1.0 fm. This can be 
satisfied e.g. by u s = 0.6 c and T = 0.12 GeV. 



C(k z ,q x ) = l + exp(-i? 2 g 2 ) x 



cosh ( ?^ 



cos (2^(4) 



cosh 



( 2~fk,v z \ 



1 



C(k z ,q y ) 
C(k z ,q z ) 



1 + cxp(-i?V) , 



:>2„2\ 



(60) 



1 + cxp(-ir^) x 
cosh(?^) + cosh(^) 



cosh 



(^) + 



1 



Therefore only the correlation functions in the k Zl q x - 
direction and in the g z -directions are affected by the z- 
directed velocity of the source. In this direction, k z , un- 
fortunately it is difficult to detect the two particle corre- 
lations. For the k x and fc y -directions the ^-distribution 



is affected by the displacement of the two sources by ±d x . 
The ^-distribution is not effected by either the displace- 
ment or the source velocities. 



2,0 
1,8 
1,6 



cr 
O 



1,4 



1,2 



1,0 



__, 1 , f=_ — , 

1 2 3 

q (1/fm) 

FIG. 6. (color online) The correlation functions, C(k,q) 
for two moving sources where the displacement of the sources 
is in the ^-direction, and the center-of-mass momentum, k, 
of emitted particles is in the x-direction. The dashed green 
lines are for the relative momentum, g z , the solid red line 
is for q y and dotted blue lines are for q x . For large values 
of the center-of-mass momentum k x the correlation functions 
C(k x , q x ) and C(k x ,q z ) will approach the correlation function 
C(k x ,q y ) (red line). For q x (blue lines) the displacements are 
d x — 1,0 fm, and for q z (green lines) the velocity is chosen 
such that 7t> z /T s = 1.0 fm. The values of k z are for the blue 
lines: 0.25, 0.5, 0.75, 1.0, 1.25 and 2.0 far 1 and for the green 
lines: 0.25, 0.5, 0.75, 1.0 and 1.5 fm" 1 . 

The correlation function for different source locations 
and velocities are similar. The cosine term appears in 
the same direction as the axis at which the sources are 
located and the hyperbolic cosine in the direction of the 
velocity. See Figs. [5]and[6j The zero points discussed for 
the two static sources at Eq. (51), appear in the distri- 
butions C(k x ,q x ) and C(k y: q x ). These distributions do 
depend on the magnitude of the flow velocity, v z , but not 
on its direction! This arises from the fact that the detec- 
tors are assumed to be reached from both sides of the 
system with opposite velocities with equal probability. 

Unfortunately the dominant direction of flow (see Fig. 
[8| is the beam direction (z— direction), where we have 
no possibility to place high acceptance detectors. At the 
same time the strongest effect of the flow appears in this 
direction. 

The rotation in the reaction plane can also be charac- 
terized with another configuration of the two moving 
sources, when the displacement is in the z-direction while 
the flow velocities are pointing into the x-direction, so 
that the source at x\ = d z has a negative velocity, — v x 
while the source at X2 — ~ d z has a positive velocity, v x . 
(See Fig. [7]) The detailed description of the correlation 
functions from this configuration can be obtained in a 
straightforward way similarly to the previous case, see 
Eq. (61 ). In this case the flow has the most dominant ef- 
fect in the A^-direction, which is accessible for detection. 
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The cc-directed flow, however, is more sensitively depen- 
dent on secondary effects, like the Kelvin-Helmholtz In- 
stability 0. 



-4 



FIG. 7. (color online) Two moving sources in the reaction 
([a; — z]) plane with a distance between them of 2d in the 
z— direction. The sources are moving in the directions indi- 
cated by the (red) arrows. 

In this configuration of the sources the magnitude of 
the flow velocity makes visible change in C(k,q), in the 
(k x , g^-direction also, which is detectable by the usual 
detector configurations. Still the direction of the rotation 
does not appear in the observables with the approach 
presented here. 

This actually arises from the simplifying assumption, 
that the freeze out is happening instantly at a timclikc 
hypersurface with ct m = (1,0,0,0), where particles from 
all sides of the system can reach each detector with the 
same probability. We will return to this problem after 
having discussed the more complex source configurations. 

C(k x , q x ) = 1 + cxp(-i? 2 ^) x 

cosh(^^f^) +cosh(») 



cosh 



C{k x ,q y ) = l + cxp(-i? 2 ^) , 
C(k x ,q z 



l + exp(-R 2 q 2 z ) 
cosh 



(61) 



cos (2q z d z ) 



cosh ( 



1 



For these two-particle correlation measurements it is 
necessary to identify independently, event by event the 
global collective reaction plane azimuth, ^rp, experi- 
mentally and the corresponding event by event center of 
mass of the system (e.g. with the method [17]). Know- 
ing these we can identify the /^-direction (and the k y - 
direction also. 

In this section we derived a relatively simple formula 
for two sources with opposite positions and opposite ve- 
locities. These kind of systems were analysed earlier for 
radially expanding systems. 

Recently due to the angular momentum in peripheral 
heavy ion collisions strong rotation [TJ and turbulence 
(Kelvin-Helmholtz Instability) [5] were predicted in fluid 



Projectile 
spectators 




Target 
spectators 



FIG. 8. (Color online) Typical orientation of the spatial 
axes in case of an ultra-relativistic heavy ion reaction shortly 
after the impact. In the configuration space the projectile and 
target appear to be flat due to the Lorentz contraction. 



dynamical models arising from the symmetries, shear and 
vorticity of the initial state. 

In the simple two source example shown in the pre- 
vious section the two sources may describe a rotation 
if the sources are at a distance from the center in the x- 
direction, X\ — (+d, 0, 0) and x-i = (—d, 0, 0), while these 
have opposite velocities pointing into the z-direction, 
tti = 7(1, 0, 0, v z ) and u 2 = 7(1, 0, 0, -v z ). 



I 



4- 
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FIG. 9. (color online) Four moving sources in the reaction 
([a; — z]) plane, one pair, si, is separated in the x— directions 
and the other, s2, is in the z— direction. The sources are 
moving in the directions indicated by the (red) arrows, ±u a i 
for the 1st pair and ±u S 2 for the other. 

It is important to mention that to detect rotation 
the accurate identification of the reaction plane and its 
proper orientation is necessary. In the so called "cumu- 
lative" methods the reaction plane is identified but its 
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projectile and target sides are not. This makes it im- 
possible to detect directed flow, and odd components of 
the global collective flow. (All harmonic components of 
random fluctuations of course can be detected.) Fur- 
thermore, not only the reaction plane with proper direc- 
tion but also the event by event center of mass (cm.) 
should also be identified [37]. This hardly ever done! In 
both cases the use of zero degree calorimeters are provide 
an adequate tool as these are sensitive to the spectator 
residues. 

The correlation function depends both on vectors k 
and q. To detect rotation the choices should be correlated 
correctly with the beam and the directed reaction plane 
as illustrated in Fig. [8] The positive x-axis points in the 
direction of the projectile, which moves in the positive 
direction along the z-axis. 



In Eq. (57 1, in the above situation, ku s 



lk z v Z: 



qu s = jq z v z and qx s = q x d. Thus, the Correlation 
function, apart of the single cell source size, R, sensitiv- 
ity, has a specific dependence on k z and q z , as well as 
on q x . Unfortunately it is difficult to measure the par- 
ticle momenta in the z-direction as it coincides with the 
beam. The q x dependence would enable us to estimate 
the distance of the two sources. 



C. Four Fluid Cell Sources 

Four sources can be treated as a combination of two 
moving double source systems. We use the same param- 
eters as under paragraph |IIB[ where s\ and S2 will be 
the two different pairs of sources with different locations 
and velocities. See Fig. [9] 

So we have 



C(k,q) = 1 +cxp(-i?V)x 



. 2k ■ u s , 
cosh 1 



2 cosh 



2 cosh 



, , q u Sl \ (2k-u S2 
cosh I — — — | cos(2q • x Sl ) + cosh ' 



k ■ (u Sl - u 
T 

k ■ (u si + u s 



. 2k ■ u s , 
cosh 1 



■- • 1 \ , / 1 ' \ U si + u s 2 ) I 

1 cosh I — I cos(q • (x Sl + x S2 ))- 



cosh 



2T S 
q ■ (u Sl - u s 



cosh 



/ 2k ■ u s 



2T, 



2 cosh 



cos(q • (tc Sl - x S2 )) 
k ■ (u Sl + u S2 



cosh | 9 ™ S2 ) cos(2q • x S2 )~ 



2 cosh 



(62) 



fc ■ (tt Sl - U s 



r 



If si = S2 then we recover Eq. (57 1 



In the case of a rotating but symmetric system the 
displacements and velocities are of equal magnitude and 
are orthogonal to each other in the two pairs: x Sl _L 
x S2 and u Sl _L u S2 . Thus a simple sign change of the 
velocity for one of the pairs or both does not change the 
result, and so the rotation can be identified, but this 
evaluation does not provide sensitivity to the direction 
of the rotation. The reason is in the simplified freeze 
out assumption as we mentioned already at the end of 
paragraph |II B| 

If the two pairs are not completely identical, i.e. the 
magnitude of the characteristic quantities of the two 
source pairs are not equal then a sensitivity to the di- 
rection of the rotation may in principle occur. However, 
if we change the direction of the velocities of the two 
source pairs simultaneously (as it happens in changing 
the direction of rotation) the result still does not change. 

Four Sources with Flow Circulation: Recent fluid 
dynamical studies indicate [U [2] , that due to the initial 
shear and angular momentum the early fluid dynamical 
development has significant flow vorticity and circulation 
on the reaction plane. These were recently evaluated [1] . 



At the present LHC Pb+Pb collision energy in the men- 
tioned fluid dynamical model calculation the maximum 
value of vorticity, w, was found exceeding 3 c/fm , and 
the circulation after 6 fm/c flow development and ex- 
pansion was still around 4-5 fm-c. This vorticity in the 
reaction plane was more than an order of magnitude big- 
ger than in the transverse plane estimated from random 
fluctuations in ref. 1281. 



In this section we will look at the four source corre- 
lation function with similar circulation as in the above 
mentioned fluid dynamical model estimates in the reac- 
tion plane. See Fig. [9| We will simulate a circulation 
value r = 5/to • c. We use Eq. ( 62 ) where the center-of- 
mass momentum, k points in the x — direction. 



Since the position and velocity arc of the same 
value and because of symmetry the correlation functions 
C{k x ,q x ) and C(k x ,q z ) provide the same values. So we 
take the correlation function C{k xi q x ) and we have af- 
terwards some simplifications. See Fig. |10| 
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FIG. 10. (color online) The correlation functions, C(k,q), 
for 4 sources where the displacement is such that there is one 
source pair on the x-axis and one the z-axis, the center-of-mass 
momentum of the emitted particles, k, is in the x-direction. 
The dotted blue lines are for the velocity v = 0.5 c and dis- 
placement d x = d z = 1.6 fm. The dashed green lines are for 
the velocity v = 0.8 c and displacement d x = d z — 1.0 fm. 
T 3 = 0.20 GeV. In both cases the circulation is T = 5 fm-c. 
The values of k x are for the blue lines: 0.5, 1.5, 3.0 and 6.0 
fm -1 and for the green lines: 0.5, 1.5 and 3.0 fm -1 . The solid 
red line is the correlation function C{k x , q y ). For large values 
of the center-of-mass momentum k x the correlation functions 
C(k x ,q x ) and C(k x ,q z ) will approach the correlation func- 
tion C(k x ,q y ). The larger displacement and smaller rotation 
velocity leads to stronger deviation from the unaffected cor- 
relation function C(k x ,q y ). 



1 + cos(2q, x d) + cosh 2 



cosh 



C{k Xl q x ) = l + exp(-.RV)x 
! ( (> h | k x yv x \ / <hTf_Vx \ cos ( ferf ) 



cosh 



2k x ^v x 



4 cosh 



kxlVx 



X 
-1 



(63) 



For C{k xi q y ) we have the same result as we had for 
the two moving sources. Here the flow and displacement 
have no effect. 

Let us look at comparisons for similar circulations and 
for similar displacements. 

By comparing Figs. [Tl] and [T2| we see that an increase 
in the displacement of the sources gives a increase in 
the apparent size of the system (narrower q— distribution. 
We also see that the measured size of the system increases 
with decreasing velocity. At the same time the shape of 
correlation functions are becoming less and less Gaus- 
sian as the flow velocities increase. At the same time the 
structure of the correlation function is also very different 
in different directions, which is not the case for spherical 
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FIG. 11. (color online) The correlation functions, C(k,q), 
for 4 sources where the displacement is such that there is 
one source pair on the x-axis and one the z-axis, the center- 
of-mass momentum of the emitted particles, k, is in the x- 
direction. The dotted blue lines are for the velocity v = 0.5 c 
and displacement d x — d z — 1.6 fm (same as in the previous 
figure). The dashed green lines are for the velocity v — 0.95c 
and displacement d x = d z = 0.84 fm, T s = 0.20 GeV. The 
values of k x are: for the blue lines 0.5, 1.5, 3.0 and 6.0 fm -1 
and for the green lines 0.1, 0.25, 0.5 and 1.5 fm -1 . The solid 
red line is the correlation function C{k x ,q y ). For large values 
of the center-of-mass momentum k x the correlation functions 
C(k x , q x ) and C(k x ,q z ) will approach the correlation function 
C(k x ,q y ). Now for the dashed green lines with even higher 
velocity and smaller displacement, the deviation is significant 
and it is in the positive direction. 



or linear expansion. This indicates that the rotating sys- 
tem contributes to essential non-Gaussian modifications, 
which can be seen directly in the correlation function, but 
they would become invisible if we would like to fit these 
data with a set of Gaussians. Earlier works studied the 
correlation function at different angles or pseudorapidi- 
ties with Gaussian parametrizations 0^], however, for 
rotating systems this is not the most sensitive way of pre- 
senting the results. Thus rotation can be detected even 
in " symmetric" few source systems where the emission is 
equally probable from all emitting sources. This emission 
scenario is less applicable to emission from heavy ion re- 
actions where the absorption of particles in QGP is not 
negligible, and this affects the emission from the interior 
of a timelike (spacelike) FO layer, where the emission of 
earlier (deeper) emitted particles are quenched. 

The correlation function is symmetric in all these cases 
as sources from opposite sides of the system contribute 
equally. Thus the correlation function is not sensitive to 
the direction of rotation. 



III. ASYMMETRIC SOURCES 

We have seen in the previous few source model ex- 
amples that a highly symmetric source may result in 
highly symmetric correlation functions, however, this re- 
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FIG. 12. (color online) Same figure as the previous one, 
however the circulation is not the same and the displacements 
are equal. The dotted blue lines are for the velocity v = 0.5 c 
and displacement d x = d z — 1.0 fm. The dashed green lines 
are for the velocity v = 0.95 c and displacement d x — d z — 1.0 
fm. T s = 0.20 GeV. The values of k x are for the blue lines: 
0.5, 1.5, 3.0 and 6.0 fm~ and for the green lines: 0.1, 0.25, 
0.5 and 1.5 fm -1 . Here the displacement is the same but the 
ultra-relativistic velocities lead still to the deviation in the 
positive direction. 



suits were not sensitive to the direction of the rotation, 
which seems to be unrealistic. We saw that this result is 
a consequence of the assumption that both of the mem- 
bers of a symmetric pair contribute equally to the cor- 
relation function even if one is at the side of the system 
facing the detector and the other is on the opposite side. 
The dense and hot nuclear matter or the Quark-gluon 
Plasma are strongly interacting, and for the most of the 
observed particle types the detection of a particle from 
the side of the system, - which is not facing the detector 
but points to the opposite direction, - is significantly less 
probable. The reason is partly in the diverging velocities 
during the expansion and partly to the lower emission 
probability from earlier (deeper) layers of the source from 
the external edge of the timelike (or spacelike) FO layer. 
This feature is recognized for a long time and discussed 
in detail by now. This influences the particle emission 
(or freeze out (FO)) process and modifies the post FO 
particle distribution. This topic has an extended litera- 
ture, and this feature destructs the symmetry of emission 
of from source pairs at the opposite sides of the system 

EDUHIinHISlEilHSE]. 

For the study of realistic systems where the emission is 
dominated by the side of the system, which is facing the 
detector, we cannot use the assumption of the symmetry 
among pairs or groups of the sources from opposite sides 
of the system. Even if the FO layer has a time-like normal 
direction, the (k^a^) factor yields a substantial emis- 
sion difference between the opposite sides of the system. 
Now we want to demonstrate this effect on few source 
examples, and we will demonstrate the consequences of 
the non-symmetric emission. 



A. The Emission Probability 

It was first recognized that the freeze out with the 
Cooper-Fry description |14j , may lead to negative contri- 
butions for particles, which move towards the center of 
the system and not in the direction out, towards the de- 
tectors. The first proposal to remedy this problem came 
from Bugaev |30j , which led to the introduction of an im- 
proved post freeze out distribution in the Cooper-Frye de- 
scription, first with the Cut-Juttner distribution 13"2] 
and then by the Cancelling- Juttner distribution [36j . 

Subsequently it was realized that for the realistic treat- 
ment of the freeze out process in transport theory one has 
to modify the Boltzmann transport equation by replac- 
ing the local molecular chaos assumption with a non-local 
one, where the point of origin is also included in the phase 
space distributions of the colliding particles. This led to 
the Modified Boltzmann Transport equation (MBT), and 
also the necessity to introduce an escape probability, P esc 
was pointed out. 

The escape probability was then introduced and anal- 
ysed in a series of publications [17U19, 38J, in transport 
theoretical approaches. It was pointed out that even if 
the pre FO distribution is a locally equilibrated isotropic 
distribution, the freeze out process and the escape proba- 
bility will provide a nonisotropic distribution which elim- 
inates the earlier observed problems. This developing 
anisotropy in the freeze out process occurs for freeze out 
both in space-like and time-like directions. 

The escape probability introduced in the works |17l - 
IT!?1 l3"8j . for a space-time surface layer of the system of 
thickness L, pointing in the four direction was given 
at a point x^ inside the freeze out layer as 

P -^ K (r^)(^) e ^> (64) 

where p M is the momentum of the escaping particle, u M (a;) 
is the local flow velocity and s — x^a^ is the distance of 
the emission point from the inside boundary of the layer. 
The first multiplicative term describes higher emission 
probability to the particles, which are emitted closer to 
the outside boundary of the layer, the second multiplica- 
tive term describes the higher emission probability for 
the particles, which move in the normal direction of the 
surface, because these should cross less material in the 
layer. The last term secures that only those particles can 
escape, which move outwards through the layer. 

The last two momentum dependent factors are impor- 
tant in transport theoretical models, to determine the 
shape of the post FO momentum distribution, e.g. |36j . 
which would replace the Juttner distribution. This shape 
modification happens to the single and two particle dis- 
tributions equally, and it acts in all emission directions, 
k, equally, so this effect is secondary from the point of 
view of the flow velocity dependence of the correlation 
function. 

In order to describe the complete freeze out process 
for a reaction the system had to be surrounded with a 
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freeze out layer in the space-time, and the phase space 
distribution of the escaping, frozen out particles can be 
obtained by integrating over the whole 4-volume of the 
freeze out layer the local (usually isotropic) phase space 
distribution with the escape probability P esc (x). This 
procedure would then play the role of function G(x) in 
the source function in Eq. ( |10[) instead of the simplified 
assumptions, as e.g. in Eq. ( |14[ ). 

The correlation function, C(k, q) is always measured in 
a given direction of the detector, fc. Obviously only those 
particles can reach the detector, which satisfy k^a^ > 0. 
Thus in the calculation of C(k, q) for a given fc- direction 
we can exclude the parts of the freeze out layer where 
fc^ < (see Eq. (10) of ref. [TO] or ref. [3D]). For 
time-like FO a simplest approximation for the emission 
possibility is P esc (x) oc k^u^x) [2"U] . 

In a model calculation we therefore have to define the 
freeze out layer also, this realistically should not include 
the whole space-time volume of the reaction. In case of 
calculating C(k, q) for a given fc we should select the rel- 
evant part of the freeze out layer, which may contribute 
to emission in the k direction. This should be a layer 
of 2-3 m.f.p facing the detector at the direction fc. This 
can eliminate the symmetric pairs of fluid cells in the 
previous calculations of the correlation function, even if 
the emission normal is timelike, because the FO particle 
from an earlier emission point in the ST has to propagate 
through the plasma for some finite time, with consider- 
able quenching. 

Therefore in the following models we should apply the 
escape probability and we should define a fc-dependent 
freeze out layer also! The most simple approximation is 
to select an emission layer from the system for a given fc- 
direction with uniform emission probability from within 
this layer. The next to most simple approximation is 
to introduce an emission probability within the layer, in- 
creasing towards the outside boundary of the layer. (Here 
it is important to mention that the spatial emission prob- 
ability should be sufficiently smooth, so that one fluid cell 
and its contribution to C(k, q), should not be effected by 
this emission probability. 

When we have up to 4 sources we can always add fc- 
dependent emission weights to these sources. This still 
would qualitatively change the outcome. As we discuss 
here up to four sources only a detailed formal evaluation 
of the emission probability would be an exaggerated ap- 
proach, by defining more parameters than the outcome, 
so we just define the weights themselves here. In a full 
3+1D fluid dynamical model with 100000+ fluid cells of 
course we have to apply a realistic and general evaluation 
of emission probability for every point of the ST. 



B. Emission probabilities for few sources 

Two sources: The previous discussion included two 
sources (i) in the beam-, z— direction and (ii) in the trans- 



verse direction in the reaction plane, a;— direction. In case 
(i) the emission could be different from the two sources 
if the detector is in the z— direction, which is difficult to 
achieve, so we do not have to discuss this possibility. 
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FIG. 13. (color online) Two moving sources in the reaction 
([a; — z\) plane, separated in the x— direction (case (ii) in the 
text). The sources are moving in the directions indicated by 
the (red) arrows. The detector is in the positive x— direction, 
thus the source on this side has more dominant emission into 
this direction, and this is indicated by the bigger size of the 
source on this side. 

In configuration (ii) the observation can be in different 
fc-directions. If fc points into the ±y— direction, then the 
probabilities must be identical so emission probabilities 
do not lead to any change. 

If fc points into the ±x— direction, then one of the 
sources is closer to the detector and may shadow the 
more distant one. Thus, we can just introduce two pos- 
itive weight factors so that w c is the weight for the cells 
closer to the detector and w s is for the cells which are 
far from the detector measuring the average momentum 
fc. These weights are the same for the calculation of the 
nominator and denominator of the correlation function, 
so their normalization does not influence the correlation 
function. 

As not all emitted particles reach a given detector the 
normalization is also dependent on the direction of the 
detector. Thus, we evaluate the correlation function this 
way. This immediately changes the earlier result ( 58 1 , be- 



cause it breaks the symmetry between the two sources. 
We can simply repeat the calculation for two moving 
sources in section II B modifying the derivation of Eq. 



( 56 ) and obtain the general result 



C(k,q) 
iku s 



1 + exp(— R q 



2ku s 



X 



2w r w* cosh 



( 93h 



cos(2qa; s 



2ku. 



2ku, 



2w r w q 



(65) 



Note that this result is valid for the case when fc points 
to the +x direction, because the weights depend on this 



and w c > w s . See Fig. 13 The fact that the emission 
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from the source, which is closer to the detector is stronger makes the direction of the flow detectable. 

If we introduce the notation w c — 1 + e and w s — 
1— e, the deviation from the symmetric result will become 
apparent 



C(k,q) 



1 + exp(--RV) 



(1 



) cosh 



2ku s 



2esinh + (1 - e 2 ) cosh (%gA cos(2qa; s ) 



+x 



(1 + e 2 ) cosh (^A + 2esinh ( 



2ku, 



+ (l-e») 



(66) 



If e — > 0, i.e. if w c = u> s , we recover the earlier result, 



Eq. (57). 
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FIG. 14. (color online) Two moving sources in the reaction 
([x — z]) plane, separated in the x— direction. The sources are 
moving in the directions indicated by the (red) arrows. The 
two "tilted" detector directions are indicated by the (blue) 
arrows labeled with fc+ and 

If e = we have the symmetric situation where both 
sources have equal contribution, the asymmetric terms 
vanish, and the result becomes to be symmetric for the 
change of the direction of the flow velocity. If e reaches 
its maximal value, e = 1 the contribution of the far side 
source is eliminated (w s = 0, w c = 2), and only the single 
nearby source contributes to the correlation function. In 
this case the asymmetric term in the nominator vanishes, 
the remaining terms in the nominator and denominator 
are equal, and we recover the single static source result. 

This result has terms, which change sign if the flow ve- 
locity, u s changes sign. The result is valid only if the de- 
tector is in the k = (1,0,0) direction. For this direction, 
however, if the flow velocity points in the z-direction, i.e. 
orthogonal to k the asymmetric term does not provide 
any contribution, so it will not show up in C(k x ,q). To 
circumvent this problem we should study detector direc- 
tions, which do not coincide with the primary axes of 
the given event (where x is the direction of the impact 
parameter vector, b, pointing to the projectile; y is the 
other transverse direction; and z is the direction of the 
projectile beam). 

Correlation in Tilted Directions: The form of the 
correlation function is the same if k is in the same plane, 




q (1/fm) 

FIG. 15. (color online) Difference of the forward and back- 
ward shifted correlation function, AC'(k±, q ut), for the value 
e = 0.50. The solid black lines are for the velocity v z = 0.5 c, 
dotted blue lines are for the velocity v z = 0.6 c and dashed red 
lines are for the velocity v z = 0.7 c. Displacement is d x = 1.0 
fm, T s = 0.139 GeV and a = b = l/y/2. The values of k are: 
for the solid black lines 0.25, 0.50, 2.00, 2.75 and 3.50 fm" 1 , 
the dotted blue lines 0.25, 0.50, 1.00, 1.75 and 2.50 fm" 1 , and 
the dashed red lines 0.25, 0.50, 0.75, 1.25 and 1.75 fm" 1 . 



the reaction plane, but it has a z component also, i.e. 
k = (k x , 0, ±fc z ). This is possible for all LHC heavy ion 
experiments, ATLAS, CMS and even ALICE, where the 
longitudinal acceptance range of the TPC (Arj < 0.8) is 
the smallest. See Fig. [Ti] 

Earlier for spherically or longitudinally expanding sys- 
tems the dependence of the correlation function on the 
tilt angle or width parameter was analysed in detail in ref. 
|10j . We do not go into similar fine details, just demon- 
strate the possibilities for an arbitrary configuration. 

Depending on the detector acceptance we should chose 
a detector direction where \k z \ is as big as the detector 
acceptance allows it. For this configuration the form of 
the correlation function is the same as ( 66 1 



C{k,q) 



C(k,q) 



+x,±z 



(67) 



+x 



with keeping the different weights, w c , w s or e so that the 
forward shifted and backward shifted directions have the 
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same weights. These weights are not specified up to now where a 2 + b 2 = 1, see Fig 
anyway. vector, q, can be measured in the directions 

For detection of the correlation function we have to 
introduce here the usual, fc-dependent coordinate system 



Then the difference 



kout = 0, 0. ± b ), 1x = a\q\, q z = ±b\q\ 
to classify the direction of q. Thus if ^side = (0, 1, 0), q y = \q\ (69) 



Qiong = (T&, 0, a), q x = Tb\q\, q z = a\q\ 
k± = (a, 0, ±6)fm _1 , k x = a\k\, k z — ±b|fe|, (68) This leads to the following correlation functions 



C(fe(±), qout) = 1 + exp(-i?V) 



(1+e 2 ) cosh (22^) +2e sinh(^fe^) + (1-e 2 ) cosh (») cos (q x d x ) 



(1+e 2 ) cosh (^^f^) +2esinh 



2-yk z v z 



+ (!-<?) 



C(k(±),q sid e) = 1 +exp(-i?V) , 



(70) 



C(k(±),qiong) = l + exp(-R 2 q 2 ) 



(1+e 2 ) cosh + 2esinh( 2 ^f^) + (1-e 2 ) cosh (W^) cos (q x d x ) 



(1+e 2 ) cosh ( 2 ^± 



2e sinh 



+ (1-e 2 ) 



Although, it seems that C(k(±^, q ut) and C(k(±^, qiong) are the same, this is in fact not the case, because the values 
of the components of the different types of k and q are not the same as described in Eqs. ( 68|69 ). In all cases, the 



out-, side- and long- q — \q\. We will also use the notation k = \k\ and jv x = u x , "fv y = u y , 7« 2 = u z , so that 
u s = (u x , u y ,u z ). For example for the out component the difference of the forward and backward shifted correlation 
functions is 



AC(k ± ,q out ) = C{k + ,q out ) - C(k-,q out ) 



4exp(— R 2 q 2 ) esinhf 



2u, bk 



(1-e 2 ) 



1 - COsh ' ^2 



cos (aqd x 



(1+e 2 ) cosh (^) + (1-e 2 ) 



4e 2 sinh 2 ( ?^ 



(71) 



As Eq. (711 and Fig. 



15 show, the Differential Cor- 
relation Function (DCF), AC(k±,q out ), is sensitive to 
the speed and direction of the rotation, and it is also 
sensitive to the amount of the tilt in the directions of 
the detection, regulated here by the parameters a and 
b. AC(k±,q out ) tends to zero both if q — > and if 
q — > co. The structure of AC(k±,q out ) is determined 
by the cosh (u z bq/T s ) cos (ad x q) product. If in both ar- 
guments the coefficients of q, u z b/T s and ad x are pos- 
itive, smaller than one, and u z b/T s < ad Xl then the 
DCF is positive. If the coefficient ad x exceeds one the 
cos function changes sign at high q out values (e.g. above 
q = 1 — 2fm~ 1 ), and the DCF becomes negative at high 
qout values. Note that the ratio of the two coefficients is 
influenced be the choice of the tilting angle, i.e. by the 
parameters a and b. 

If the parameter ad x remains constant, about 1 fm, and 
then when u z b/T s becomes larger (than one) the Differ- 
ential Correlation Function becomes negative at small q 
values. 

If the parameter u z b/T s remains constant, and about 
1 fm, then when ad x becomes less (than one) the DCF 
becomes negative at small q values. 

In case if the detector has a narrow pseudorapidity 
acceptance, then k± is close to k x , i.e. b <C a and then the 
weights are maximal for the source in the x— direction, 
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FIG. 16. (color online) Four moving sources in the reaction 
([x — z]) plane, separaed in the x— and z— directions. The 
sources are moving in the directions indicated by the (red) 
arrows. The "tilted" detector directions are indicated by the 
(blue) arrows. 



as indicated in Fig. [TBJ 

If we change the direction of rotation to the opposite 
the Differential Correlation Function changes sign due to 
the sinh function in the nominator. In this configuration 
with the change of the tilt of the detector directions we 
can adjust the DCF, to the threshold value where the 
AC(k±,q out ), is still positive, which provides a sensitive 
estimate for the rotation velocity at Freeze Out. 
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This very sensitive behaviour is rather special and it 
appears in this special two source model this way. With 
an increased resolution and with more source elements 
this strong and specific structure will be smoothed out 
to some extent. 



0,08 
0,04 



5 0,00 
O 

< 

-0,04 



-0,08 



/ 


\ 


/ 

- 'X 


\ 

a \ 
\ x 


i //,--• 


^\V\a\ . .-■ 












1 1 1 1 



1 2 3 

q (1/fm) 

FIG. 17. (color online) Correlation function difference for 
the weight factors: uj a — 1.25, ut = 0.75, uj c — uud = 1-00 
for sources placed at +x, —x, +z and —z respectively. This 
weight distribution corresponds to the configuration shown 
in Fig. |16| The solid black lines are for the velocity v z = 
0.5 c, and the dashed red lines are for the velocity v z — 0.7 c. 
Displacement is d x — d z = 1.0 fm, T s = 0.139 GeV and 
a = b — l/\/2. The values of k are: for the solid black lines 
0.10, 0.50, 1.00, and 2.00 fm -1 and for the dashed red lines 
0.10, 0.50, 1.00, and 2.00 fm" 1 . The difference is larger for 
smaller values of k. 

The sinh(2w 2 bk/T s ) term changes sign in the nomina- 
tor when u z changes sign the difference of the two corre- 
lation functions, AC(k±,q out ) changes sign also because 
all other terms are symmetric to the sign change of the 
velocity. 

This is an important observation as we can detect the 
direction and magnitude of the rotation in the reaction 
plane. This difference is also increasing with the longi- 
tudinal shift, b, of the average momentum vector, k, so 
that detectors with larger pseudorapidity acceptance can 
detect the rotation better. 

In order to perform this measurement, one has to de- 
termine the global reaction plane (e.g. from spectator 
residues in the ZDCs), and determine the projectile side 
of this plane. Furthermore the event by event center of 
mass should also be identified (using e.g. the method 
shown in ref. |27j). This will be the positive x-direction. 
Then the correlation function can be measured for four 
different fc-directions in the global reaction plane. These 
four directions are shifted forward and backward from 
the center of mass symmetrically on the projectile side, 
and there should be a symmetric pair of detection points 
in the target side of the reaction plane too. 

The k directions opposite to each other across the 
cm. point give the same result, while the difference, 
AC(k±,q out ), between the Forward (F) and Backward 



(B) shifted contributions will characterize the speed and 
direction of the rotation. This symmetry can be used to 
eliminate the contribution from eventual random fluctu- 
ations. The observed F/B asymmetry depends on the 
parameters e, v z and d x , these can be estimated by mea- 
suring the correlation functions at all possible moments 
k. 

Fig. [15] indicates that the differential correlation func- 
tion has a larger amplitude for smaller k values, and the 
zero points are sensitively dependent on the rotation ve- 
locity. 

The zero points come from the term 



1 — cosh 



u z bq 



cos(aqd x ) = 



(72) 
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FIG. 18. (color online) Correlation function difference for 
the weight factors: uj a — 1.25, cut = 0.75, uj c — cod = 1-00 
for sources placed at +x, —x, +z and — z respectively. This 
weight distribution corresponds to the configuration shown 
in Fig. |16| The solid black lines are for the velocity v z = 
0.5 c, and the dashed red lines are for the velocity v z — 0.7 c. 
Displacement is d x — d z = 1.0 fm, T s = 0.139 GeV and 
a = 4b — 4/y/Vj. The values of k are: for the solid black lines 
0.10, 1.00, 3.00, and 5.00 fm -1 and for the dashed red lines 
0.10, 1.00, 3.00, and 5.00 fm" 1 . For a detector with narrow 
pseudorapidity acceptance we could use a value a > b. For 
this figu re we have used a = 46. We can compare this with 
Fig. 17 and we see that for q > 1.5 fm -1 they are very similar. 
The DCF has a smaller amplitude in the positive direction 
for values around q < 1.0 fm -1 , but a larger amplitude in 
the negative direction for values around 0.75 fm -1 < q < 1.5 
fm - . We also see that the DCF is not close to zero for values 
around 2.0 fm" 1 < k < 5.0 fm" 1 . 

a nd it is not dependent on fc, so for the values used in 
Fig. 15 and q — xfm" 1 we have 



cosh 



0.197 v z /c x 
0.139^1^17^72" 



(73) 



Since the cosine term must be positive, there are no zero 
points for tt/v^ < x < 3tt/\/2 or 2.23 < x < 6.66. 
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For x > 37r/\/2 the exp(— R 2 x 2 ) term will be small and 
there would be no correlation difference. So we will look 
at the values < x < ir/V2 

If v z is smaller than 0.57c then there will be no zero 
point for < x < tt/v2- For v z larger than 0.58c there 
will be 1 zero point for < x < ir/V2. 

The zero point for a given velocity can be found by 

0.6c 
and 



solving equation ( 73 ) numerically. For velocities v 
and v = 0.7c we have zero points at q = 0.83 /m 
q = 1.69 fm respectively. 

This indicates the sensitivity of the method and the 
possibility to influence it by the choice of the detector 



directions (via the choice of a and b) . 



C. Emission from four sources 

With four sources we can illustrate the possibilities of 
differential HBT method studies in different directions. 
The correlation functions can be calculated in general for 
four sources and two detector positions. This can then 
be applied to different detector configurations. 

The out component of the four source correlation func- 
tion with weight factors u ai ui b , u) c , uj d is given by 
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(74) 



Two examples on different detector configurations are 
given in Figs. [16] and [19] We use the same equations as 
in the two source model, Eqs. (68) and (69). 



A source with a larger weight factor is closer to the 
detector, so that u a , oj b , u> c , uid correspond to x s = 
(r x ,r z ) = (d x ,0),(-d x ,0),(0,d z ),(0,-d z ) respectively. 

In case if the detector has a wide pseudorapidity ac- 
ceptance, then k± can deviate significantly from k Xj i.e. 
b > a and then the weights are maximal for the two 
sources closest to fc + or k as indicated in Fig. 19 

Eq. ( 74 ) can be used to find the difference of the for- 
ward and backward shifted correlation function. We will 
use that d x = d z , v x = v z and a = b. 

Some examples for the differential correlation functions 
are shown in Figs. 17 and 20 Here due to the simpli- 



fied few source model we specified the weight distribution 
among the sources in a simplified way. For realistic high 
resolution fluid dynamical model calculatios the realistic 
evaluation of emission probabilities is necessary. 

We can compare Fig. [T7] with the previously shown 
two source model, Fig. |15| and we see that the ampli- 
tudes are similar but the shapes are different. First of 



all the sensitivity on the direction of rotation remained 
the same as in the simpler two source model. The two 
extra sources, c and d lead to higher amplitude for the 
Differential Correlation Function, while the regular posi- 
tions of the locations of the zero points are varying due 
to more sources with different weight parameters. 

In Fig. 18 we show the DCF for a configuration where 
the deviation between k + and fc_ is smaller, like shown 
in Fig. |16| This configuration can be applied in detectors 
where the pseudorapidity acceptance range of the detec- 
tor is not wide. Still the rotation is well detectable. In 
this configuration the accurate determination of the reac- 
tion plane and the participant center of mass momentum 
is more important. 

If the detector acceptance is wider, then the two detec- 
tors can be placed at more different angles. This config- 
uration makes the forward and backward placed sources 
more accessible to the forward and backward detectors, 
respectively. This is taken into account in the emission 
weights of our sources. These weights are now different 
for the two components of the DCF! 

The result shows the tendency that the DCF has a 
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FIG. 19. (color online) Four moving sources in the reaction 
([a; — z]) plane, separated in the x— and z— directions. The 
sources are moving in the directions indicated by the (red) 
arrows. The "tilted" detector directions are indicated by the 
(blue) arrows. In the two configurations, (a) and (b) the 
detector directions are different and the weights of the sources 
are also different, so that the sources closer to the detector 
direction have larger weights. 



similar structure in the two source model and the four 
source model in a resembling configuration. Fig. 20 has 
the same shape as Fig. fT7l but the amplitude is larger. 



For a set of large number of sources, forming a system 
with close to perfect rotational symmetry, a single cor- 
relation function would not depend on the (polar) angle 
of the detection, and the DCF would vanish. Thus, the 
DHBT method would not be applicable for highly sym- 
metric systems, like for a rotating star observed from 
within the plane of the rotation. At the same time for 
a rotating binary star system the DHBT method would 
work. Also the weighting of the sources should be differ- 
ent: If the observer is in the plane of rotation, the distant 
star is shadowed by the front one at some periods, just 
like emission from a highly opaque plasma (evidenced 
by jet quenching). If the observer is slightly out of the 
plane of rotation then the two stars are visible all the 
time and then the (time dependent) correlation function 
would change between the configurations of Figs. 4 and 
7. This also illustrates the role of symmetric and asym- 
metric weightings. 

The rotating and expanding final state of a relativistic 
heavy ion reaction is of course does not look like a per- 
fect wheel, so the four source model is a more adequate 
approximation than a wheel would be. 



IV. CONCLUSIONS 

In this work we attempted to study the possibility of 
detecting and evaluation the rotation of a source by the 
specific use of the Hanbury-Brown and Twiss method for 
rotating systems. Our primary interest was the appli- 
cation for ultra-relativistic heavy ins where in peripheral 
collisions at ultra-relativistic energies the system can gain 
large angular momentum. Nevertheless, some of the con- 
clusions can be applied to macroscopic systems also, like 
for past rotating stars. 

We selected one of the several methods to evaluate two 
particle correlations, which was suitable to study collec- 
tive fluid systems with significant and well defined inter- 
nal fluid dynamical motion. The obtained standard cor- 
relation functions were showing the consequences of the 
flow, but for highly symmetric sources the correlation 
functions gave symmetric results, which were invariant 
for the change of the direction of rotation. 

It turned out that it is important to take into account 
that the particles reaching the detector cannot reach it 
with equal probability from the near side and the far side 
of the emitting object. With this fact considered we could 
obtain correlation functions, which reflected the proper- 
ties and also the direction of the flow. These results can 
be used rather generally. 

The obtained results have shown that the correlation 
function is most sensitive to the rotation if it is measured 
in the beam direction (or close to it). This, unfortunately, 
is not possible in most heavy ion accelerator experiments, 
so we introduced and investigated a Differential Hanbury 
Brown and Twiss method, which made it possible to trace 




FIG. 20. (color online) Correlation function difference for the 
weight factors: ui a = 1.25, Ub = 0.75, uj c = 0.50, tod = 1-50 for 
sources placed at +x, —x, +z and — z respectively (Fig. |19fc ) 
and uj a = 1.25, u)t = 0.75, oj c = 1.50, U)d = 0.50 for sources 
placed at +x, —x, +z and —z respectively (Fig. 19 j). The 
solid black lines are for the velocity v z = 0.5 c, and the dashed 
red lines are for the velocity v z = 0.7 c. Displacements are 
d x = d z = 1.0 fm, T s = 0.139 GeV and a = b = 1/V5. The 
values of k for both series of lines are 0.10, 0.50, 1.00, and 
2.00 fm . The difference is larger for smaller values of k. 
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down the rotation in relativistic heavy ion collisions by 
measuring the correlation functions in the reaction plane 
at nearly transverse angles to the beam direction. The 
method is promising and can be performed in most heavy 
ion experiments without difficulties, as well as it can be 
implemented in different reaction models, like fluid dy- 
namical models, microscopic transport models and hy- 
brid models. In full scale theoretical models, the emis- 
sion probabilities from the FO layer have to be consid- 
ered. From the general formulas derived in the beginning 
of the paper apparently these dependencies can be fac- 
torized. 

To complement these studies we also applied the 
method to a high resolution, 3+ ID, computational fluid 



dynamics model [35], which was used earlier to predict 
roation, KHI, flow vorticity, and polarization [T| 12) 14) 140] . 
The result shows that the method can detect rotation, 
while the effects of irregular shape, sperical flow, and 
specific flow patterns, require a more extended analysis 
to separate all these effects form one another. 
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